hep-lat/9809045 
September, 1998 



1 



00 
On 

o\ 

<D 

in 

(N 
(N 

> 

in 
o 

On 
O 
0O 
On 



X 



What's new with the electroweak phase transition?* 
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We review the status of non-perturbative lattice studies of the electroweak phase transition. In the Standard 
Model, the complete phase diagram has been reliably determined, and the conclusion is that there is no phase 
transition at all for the experimentally allowed Higgs masses. In the Minimal Supersymmetric Standard Model 
(MSSM), in contrast, there can be a strong first order transition allowing for baryogenesis. Finally, we point out 
possibilities for future simulations, such as the problem of CP-violation at the MSSM electroweak phase boundary. 



1. INTRODUCTION 

Primordial nucleosynthesis computations tell 
that the net baryon to photon number ratio r] 
in the Early Universe is a definite non-vanishing 
number, r\ — (1...9) x 1CP 10 . On the other hand, 
it is natural to assume that after inflation, 77 = 0. 
The last instance in the post-inflationary history 
of the Universe during which 77 > could have 
been generated, is the electroweak phase transi- 
tion, at T c ~ 100 GeV §. Thus, electroweak 
baryogenesis is in a way the most conservative 
scenario of baryon number generation, and at the 
same time, the only scenario which is experimen- 
tally testable in existing collider experiments. 

The scenario of electroweak baryogenesis has 
quite a few different ingredients. To generate 
a baryon number, one needs anomalous baryon 
number violating processes in the symmetric 
high temperature phase, microscopic C- and CP- 
violation, and a thermal non-equilibrium (for a 
review, see Q). It is perhaps surprising that 
many of these ingredients can be studied non- 
perturbatively with lattice simulations. Indeed, 
baryon number violation has been studied both 
in the symmetric and broken phases of the the- 
ory and references therein). Something can 
perhaps also be said about CP-violation in the 
MSSM (Sec. |). Finally, whether there is non- 
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equilibrium or not, depends on the order and 
strength of the transition. 

The main interest here will be on the last of 
these questions. In addition to what kind of 
non-equilibrium phenomena there can be, this 
determines whether the baryon number violat- 
ing processes are sufficiently switched off after 
the transition for any baryon number possibly 
produced to remain there. Indeed, while the 
baryon number violating rate is assumed to be 
only parametrically suppressed before the tran- 
sition, B ~ a^y \n(C / aw)^ (5), it is expo- 
nentially suppressed after the transition, B ~ 
exp[-45(7j ff /T)]T 4 §. The general constraint for 
the baryon number generated to remain there af- 
ter the transition is vh /T c > 1 fl],|| , and we would 
hence like to compute this ratio. 

2. PERTURBATION THEORY 

In principle, vh/T can be computed in per- 
turbation theory. So why does one need simula- 
tions? This may seem like quite a relevant ques- 
tion, as we are studying the electroweak sector 
of the Standard Model, for which perturbation 
theory works perfectly at zero temperature. 

However, things are different at finite temper- 
atures. This is due to the so-called infrared prob- 
lem. The point is simply that there is a new scale 
T, and there can thus be new expansion parame- 
ters such as g 2 T / '(Airm) , where g 2 represents the 
couplings of the theory and m the masses. In 
practice, this kind of an expansion parameter can 
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emerge from the Bose-Einstein distribution ni,: 

^2 



g nb(m) 



9" ™«t g^T_ 

e m/T _l " m 



(1) 



Thus, it is the bosonic degrees of freedom which 
are particularly problematic. Now, if m repre- 
sents, e.g., the W mass m\y — gvaj^, then m w 
in the symmetric phase, and perturbation theory 
need not work at all. The phase transition takes 
place between the symmetric and broken phases, 
and thus its characteristics may also be unreliably 
described. Thus we need lattice simulations. 

3. NON-PERTURBATIVE METHODS 

3.1. 4d simulations 

In principal, the most straightforward way 
to attack the problem is to do standard four- 
dimensional (4d) finite temperature lattice sim- 
ulations. However, in the present context 4d sim- 
ulations turn out to be quite demanding. This 
is, somewhat surprisingly, due to the fact that 
the coupling is weak. A weak coupling makes the 
system have multiple scales, and if a lattice with 
spacing a and extent N is to describe the infinite 
volume and continuum limits, one must require 
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For small g, lattices thus need to be very large. 

It is quite remarkable that in spite of this se- 
vere requirement, a continuum extrapolation can 
sometimes be carried out (see below). In partic- 
ular, one can employ an asymmetric lattice spac- 
ing H , which should essentially relax the leftmost 
inequality in Eq. 

The second problem with the 4d simulations 
is that only the bosonic sector of the Standard 
Model can be studied, since chiral quarks (espe- 
cially the top) cannot be put on the lattice. 

3.2. 3d simulations 

Another possibility is simulations in a three- 
dimensional (3d) effective theory. The main idea 
of this approach is to combine the best parts of 
perturbation theory and simulations: one can in- 
tegrate out massive modes pcrturbatively, which 
works well since the couplings are small, and then 



study light modes non-perturbatively. In the first 
step, the original 4d theory reduces to a 3d one. 

The 3d approach allows to overcome the two 
problems of 4d simulations mentioned above. In- 
deed, since the large mass scales (ttT, \f2gT) are 
removed, it is easier to satisfy Eq. (^). Conse- 
quently, an infinite volume and continuum ex- 
trapolation with a relative error of, say, 5%, can 
be obtained with rather moderate computer re- 
sources (for a review, see M). Moreover, one 
can study the full theory with chiral fermions ^ , 
for realistic (small) gauge couplings. In fact, the 
same simulation results tell about the infrared 
properties of many other theories, as well, such 
as the MSSM in a part of its parameter space |J. 

The main question concerning dimensional re- 
duction is, of course, how accurate the effec- 
tive theory constructed really is. Parametrically, 
the error 5G for static Green's functions G is 
SG/G <0(g 3 ) ||, but an essential point is to 
convert this parametric estimate into a numeri- 
cal one. What one can do is to compute some 
particular higher dimensional operators, and see 
what kind of an effect they give in different in- 
frared Green's functions. This leads to errors on 
the ~ 1% level, and often even smaller in theories 
without the top quark ||. A conservative error 
estimate is then < 5% in the Standard Model (for 
m H < 250 GeV). In the MSSM the errors can be a 
bit larger, since the strongly interacting squarks 
may play a significant role, and since there are 
more mass parameters which may compromise 
the high temperature expansion. It is perhaps 
worth stressing that perturbative dimensional re- 
duction does not work at all for the QCD phase 
transition, where the coupling is large. 



4. SIMULATION RESULTS 

We now review the main recent physics results 
obtained with the two approaches. 

4.1. The Standard Model 

In the 3d approach, the Lagrangian relevant for 
the Standard Model is 

£ 3 d = jF-F- + (D l cj ) yD^+m 2 3 ^ ( f>+X 3 (^ ( t>) 2 .(3) 
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Figure 1. The phase diagram of the Standard 
Model. The non-perturbative endpoint location 
has been studied with 3d simulations in [11-14] 
and with 4d simulations in [15-18]. In perturba- 
tion theory (dotted line), the transition is always 
of the first order. 



The U(l) group has here been neglected (i.e., 
sin 2 9w — 0), since its effects are small 0. Let 
us denote 

x = ^/gl y = rnl{gl)/gt. (4) 

In the 4d simulations, one studies the 
SU(2)+Higgs theory, whose Lagrangian is pre- 
cisely Eq. (||) but in 4d. 

The theory in Eq. (J3j> has a first order phase 
transition for small Higgs masses (small values 
of a;) 0. The transition gets weaker for larger 
Higgs masses, and ends at ran ~ 80 GeV Jll[, see 
Fig. [|. Recently, the interest has been in studying 
the endpoint region in some detail. Here, pertur- 
bation theory does not work at all and the dy- 
namics is completely non-perturbative. 

The fact that there is an endpoint, was first 
reliably demonstrated in The endpoint 

location was determined more precisely in Jl3[ . 
A continuum extrapolation of the endpoint loca- 
tion was made in E3], employing improvement 



formulas derived in IT3] : 

x c = 0.0983(15), y c = -0.0175(13). (5) 

In p^ ], it was also shown that the endpoint be- 
longs to the 3d Ising universality class. 

The values in Eq. (||) can be converted to the 
endpoint locations in different 4d physical the- 
ories, using the relations derived in ||. Some 
values are given in Table 1. The errors here rep- 
resent the errors in Eq. (^): no additional errors 
have been added from dimensional reduction. 

With 4d simulations, the endpoint location in 
the SU(2)+Higgs model has been studied at a 
fixed (symmetric) lattice spacing in [[l5|Jl6| , and 
with an asymmetric lattice spacing in |l7| , |l8[ . 
A continuum extrapolation has been carried out 
in |fl8|| , and that result is shown in Table 1. It 
should be noted that the exact MS gauge cou- 
pling to which the 4d simulations correspond, is 
not known. This affects strongly the critical tem- 
perature (T c oc mn/g), while the endpoint loca- 
tion itself is not that sensitive. 

We can now compare the 3d and 4d results for 
SU(2)+Higgs. Clearly, they are completely com- 
patible. 

Finally, consider the effect of sin 2 9w- In 
general, the hypercharge U(l) group makes the 
transition slightly stronger, though not by very 
much pc[ |. Thus one might also expect that the 
endpoint location changes to somewhat larger x 
than in Eq. (||). The infinite volume and contin- 
uum extrapolation of the endpoint location has 
not been determined with sin 2 9w = 0.23, but 
finite volumes have been studied in pOfl , On a 
lattice with A/(g 2 a) = 8 and volume = 32 3 , we 
get 

x° c = 0.1043(22), y° c = -0.02860(99), 

x\ = 0.1045(14), y\ = -0.02125(76), (6) 

where (°) refers to sin 2 9w — and (*) to 
sin 2 6w = 0.23. Hence x c does not appear to 
depend significantly on sin 2 9\y , while y c changes 
a bit. Assuming that the same pattern remains 
there at the infinite volume and continuum limits, 
the endpoint location in physical units is given in 
Table 1 also for sin 2 9 W = 0.23. 

Recent topics of interest, other than the end- 
point location, include the excitation spectrum 
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Table 1 

The endpoint location in different theories. No dimensional reduction error estimates have been added to 
the 3d results here. The values cited for sin 2 9w — 0.23 assume that the finite volume and lattice spacing 
effects do not depend on sin 2 9 Wl see the text. The comparison of critical temperatures between 4d and 
3d is sensitive to the relation of g- 2 -^ and g^, which is not known, while to# c is much less sensitive. 



Method Theory 




Couplings 


Endpoint location 




sin 2 6\v 




m H , c /GeV 


T c (m H . c )/GeV 


3d Standard Model 





0.426 


72.2(7) 


110.1(8) 


Standard Model 


0.23 


0.426 


72.3(7) 


109.2(8) 


SU(2)+Higgs 





0.57(2) 


65.9(7) 


133(4) 


4d || SU(2)+Higgs 





Sk * 0.57(2) 


66.5(14) 


128(6) 



around the endpoint |2l| , and the behaviour of 
some topology related observablcs [ p2] |. 

All in all, we can summarize the main lattice 
results for the Standard Model as follows: 
Perturbation theory: The non-perturbative 
transition is weaker than in perturbation theory 
and ends at mjj = 72(2) GeV. 
Cosmology: Although in principle all the ingre- 
dients for baryogenesis are there already in the 
Standard Model, in practice one needs something 
more, to have a first order transition. 
Dimensional reduction: In the first order 
regime, the results from 4d and 3d have been ob- 
served to be completely compatible Now a 
similar agreement has been demonstrated in the 
non-perturbative endpoint regime Q. Thus we 
can be confident that the accuracy estimates of 
dimensional reduction are reliable. This is im- 
portant since theories with chiral fermions can 
presently only be studied with the 3d approach. 

4.2. MSSM 

In contrast to the Standard Model, baryogene- 
sis could in principle work in the MSSM. This is 
because (1) there can be more CP- violation in the 
MSSM, due to new complex phases in the scalar 
sector (see Sec. ||), and because (2) there can be a 
stronger transition in the MSSM, due to a larger 
number of light bosonic degrees of freedom, viz. 
the squarks. Thus MSSM is a natural candidate 
for electroweak baryogenesis. 

The parameter regime relevant for MSSM 
baryogenesis has been studied extensively in per- 
turbation theory. The relevant case has been 
found to be that the right-handed stops are light, 



mi R < m top , and the Higgs mass is anything al- 
lowed by experiment and by the MSSM, to# = 
90. ..110 GeV [24-27]. 

As was the case for the Standard Model, per- 
turbative estimates are nevertheless not necessar- 
ily reliable. Thus, the transition has been studied 
with 3d simulations |28| ]. 

The light degrees of freedom appearing in the 
3d effective theory are now the spatial compo- 
nents of the SU(2) and SU(3) gauge fields, the 
right-handed stop U, and the combination H of 
the two Higgses Hi , H2 which is light at the phase 
transition point. The corresponding action is 

£ 3d = If*!* + \( :,)(;,) + l3 H^HU^U (7) 

+{D?H)\D?H) + m 2 H3 H^H + X H3 (H^ H? 
+{D s i U)\D s i U) + m 2 U3 U^U + X U3 (U^U) 2 . 

Here Df ,Df are the SU(2) and SU(3) covariant 
derivatives, and F°-,G^ are the corresponding 
field strengths. The parameters of this action can 
be determined in the standard way ^sjj. 

To simulate the theory in Eq. ((7J), is straight- 
forward but technically demanding. For details, 
we refer to p8| . 

The basic result of the simulations is shown 
in Fig. 0. There we show the Higgs field expec- 
tation value after the transition, for vtlh = 95 
GeV, as a function of a parameter % which de- 
termines the zero temperature right-handed stop 
mass through ~ (m 2 op — to 2 /) 1 / 2 w 160. ..150 
GeV. The lattice results are compared with 2-loop 
perturbation theory (the solid line). 

Based on Fig. H, we can summarize the lattice 
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m Q = 1 TeV, tanp = 3 {m„ = 95 GeV) 




1.0 

50 55 60 65 70 75 

Figure 2. The Higgs field expectation value in the 
MSSM f28fl , compared with 2-loop perturbation 
theory. For rhjj > 64. ..67 GeV, there is a two- 
stage transition, and the values of vh/T shown 
represent the situation after the latter stage. 



results for the MSSM electroweak phase transi- 
tion as follows: 

Perturbation theory: In contrast to the Stan- 
dard Model, the electroweak phase transition in 
the MSSM can be stronger than in 2-loop pertur- 
bation theory, even for large mn (although the 
difference is not very large). Thus perturbation 
theory gives a conservative estimate. 
Cosmology: If the non-perturbative strength- 
ening effect remains there also for larger Higgs 
masses than shown in Fig. |2|, then Higgs masses 
up to 105. ..110 GeV are allowed for baryogenesis, 
provided that m jR < m top |25|,|26| . 

5. FUTURE PROSPECTS 

As we have seen, lattice simulations have been 
quite successful in solving the problem of the elec- 
troweak phase transition: the case of the Stan- 
dard Model is now completely understood, and 
similar techniques have been applied also in the 
MSSM. One has also been able to demonstrate 



with 4d simulations that the dimensionally re- 
duced 3d theory works for non-perturbative quan- 
tities as well as for perturbative ones: the non- 
perturbative infrared dynamics is really three- 
dimensional, as it has to be (see also p9|). 

What is it then that still remains to be done? 
We list here a few open questions, related in par- 
ticular to the MSSM: 

1. What happens for other parameter values in 
the MSSM (a larger ran, non- vanishing squark 
mixing, etc), when the transition gets weaker? 

2. Could 4d simulations be used to estimate 
non-perturbatively the accuracy of dimensional 
reduction in a theory similar to the MSSM? Could 
4d simulations be used in some regions of the pa- 
rameter space ffflf where the theory is not weakly 
coupled at zero temperature, and thus dimen- 
sional reduction does not work? 

3. Finally, can one say something about CP- 
violation with simulations? The existence of CP- 
violation is one of the ingredients for electroweak 
baryogenesis: otherwise, one produces the same 
amounts of baryons and anti-baryons, and no net 
asymmetry arises. Let us discuss the last question 
in some more detail. 

To get enough CP-violation for producing the 
observed baryon asymmetry, turns out to be a 
non- trivial requirement It is hence an inter- 
esting prospect that there are new sources of CP- 
violation in the MSSM, related to the trilinear 
couplings of Higgses and the squarks. 

An intriguing observation is now that since 
one has two Higgses in the MSSM, the effect of 
the new CP-violating parameters can also prop- 
agate to a CP-violating phase between the Hig- 
gses. Thus, the idea arises that maybe the there 
is a non-trivial profile of the CP-violating phase 
at the phase boundary between the broken and 
symmetric phases [31-33]. The phase boundary 
is the region relevant for baryogenesis, and this 
scenario might allow for sufficient CP-violation, 
without violating the constraints that have to be 
satisfied in the broken phase. 

It appears that whether such a phenomenon 
takes place, can again be studied non-perturba- 
tively with a simple 3d theory. The theory in 
Eq. ([?]) is not enough, though: in that case, 
we were interested in the strength of the tran- 
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sition, and it was sufficient to study a particu- 
lar linear combination of Hi , H 2 . Now we are 



interested in a question for which both of the 



Higgses Hi , Hi should be kept in the effective 
theory. The effective theory will then be more 
complicated than in Eq. (]?]), but can be derived 
with precisely the same methods. In such a the- 
ory there are new CP-violating operators, e.g. 
1uyH\H2, which could have a non-trivial profile 
at the phase boundary between the symmetric 
and broken phases. This problem could, in prin- 
ciple, be solved with lattice simulations. 
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